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Abstract 

Qj . A nonlocal gravity model which does not assume the existence of a new dimensional pa- 

j-^l rameter in the action and includes a function with □ the d'Alembertian operator, is 

studied. By specifying an exponential form for the function / and including a matter sector with 
a constant equation of state parameter, all available power-law solutions in the Jordan frame 
I are obtained. New power-law solutions in the Einstein frame are also probed. Furthermore, the 

Cf^ ' relationship between power-law solutions in both frames, established through conformal trans- 

■ formation, is substantially clarified. The correspondence between power-law solutions in these 

two frames is proven to be a very useful tool in order to obtain new solutions in the Einstein 



frame. 

1 Introduction 

The acceleration of the Universe expansion is presently supported by a large number of independent 
sets of observational data, of very different kind [H El El HI El El [7] . Modern cosmological surveys allow 
astronomers to obtain increasingly accurate joint constraints on the set of cosmological parameters 
(see, e.g., The usual assumption that General Relativity (GR) is the correct theory of gravity 
at all scales leads to the remarkable conclusion that about seventy per cent of the energy density of 
the Universe at present must be smoothly distributed under the form of a slowly varying cosmic fluid 
with negative pressure, called dark energy. Ordinarily, in order to specify a would be component 
of the cosmic fluid use is made of its equation of state (EoS), namely a phenomenological relation 
between the pressure, p, and the energy density, p, corresponding to the considered component, e.g. 
p = WmP, where Wm is the EoS parameter. Contemporary experiments provide strong support to 
the statement that the dark energy EoS parameter is presently very close to —1. If this number 
were the exact value, this would lead us back to GR with a cosmological constant (and nothing else). 
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but a small deviation from this conclusion cannot be excluded by the most accurate astronomical 
data available. Moreover, the sign (positive or negative) or the tendency (e.g., the derivative) of this 
deviation cannot be clearly determined at present. This makes room for a number of theoretical 
models, derived from quite different fundamental theories, which can accommodate such situation. 

Actually, a few types of models exist which are able to reproduce the observed late-time cos- 
mic acceleration. The simplest, and most popular of them is ACDM, which fits a wide range of 
cosmological data [5]. In this model the dark energy component is just the cosmological constant 
which is added to the action corresponding to GR. Other models introduce a dynamical dark energy 
characterized by a varying EoS dark energy parameter, the standard way to obtain an evolving EoS 
parameter being the addition of scalar fields to the cosmological model. Actually, the evolution of 
the Universe is sufficiently well described by cosmological models with scalar fields, in particular, 
by quintom models, which involve two of them: a phantom scalar field and an ordinary scalar (see 
e.g. [9]). Quintom models are being very actively studied at present [TOl [TTl [121 1131 E] (Jot reviews, 
see also [15]). We should note, however, that the origin of this fluid, which produces anti-gravitational 
effects, still remains a mystery. Other popular theories involve modifications of Einsteinian gravity, 
as for instance F{R) gravity, with F{R) an (in principle) arbitrary function of the scalar curvature 
R (for recent reviews see, e.g., [161 HH IIHI [El ISUl [2ll |22]). 

Higher-derivative corrections to the Einstein-Hilbert action are being actively studied in the 
context of quantum gravity (as one of the first papers on this subject we can mention p3j; see also [21] 
and references therein). A nonlocal gravity theory obtained by taking into account quantum effects 
was proposed in [2S]. Also, string/M theory is usually considered as a possible frame (expectedly, 
the ultimate one) for the discussion of all fundamental interactions, including gravity, consequently, 
the natural appearance of non-locality within string field theory provides a very strong motivation 
for studying nonlocal cosmological models. It should be emphasized in this context that most of the 
nonlocal cosmological models available explicitly include a function of the d'Alembertian operator, □, 
and either directly define a nonlocal modified gravity I251[271[2Sl[2a[3ni[SIl[221[331[3a[33[2ni[23 
or, alternatively, add a nonlocal scalar field, minimally coupled to gravity [39] . 

In the present paper, we consider a nonlocal gravity model which contains a function of the 
operator but does not assume the existence of a new dimensional parameter in the action [25j. For 
this kind of nonlocal models, an explicit technique for choosing the distortion function, /(D"^), so 
as to fit an arbitrary expansion history, has been derived in [21], and the specific nonlocal model 
considered has a local scalar-tensor formulation [27]. The perturbation analysis of this model has 
been carried out and the Solar System test has been performed in [30]. De Sitter solutions and 
expanding universe solutions with a ~ have been investigated in [271 EH [351 ES] . In [29] the ensuing 
cosmology describing the four basic epochs was studied for nonlocal models involving, in particular, an 
exponential form for the function f{ri). An explicit mechanism to screen the cosmological constant 
in nonlocal gravity was discussed in [321 [Ml [3S]- framework of the local formulation, a 

reconstruction procedure was proposed in [29] and it has been developed in [361 [371 [38] . 

The example most usually studied [271 [2H1 [291 [331 [311 [351 ES] of a model of this kind is characterized 
by an exponential function f(n~^R) = foe°'^° ^\ where /o and a are real parameters. In this paper, 
we explicitly consider the case of an exponential function /. 

Conformal (Weyl) transformations are widely used in scalar-tensor theories of gravity, the theory 
of a scalar field coupled nonminimally to the Ricci curvature, R, and in F{R) gravity theories [IHl 
[HI [12] (see also [TH [20l [21])- The Hilbert-Einstein action and the modified gravity action can be 
related by conformal transformation [TH [201 [211 [121 [HI [B], being the corresponding equations also 
connected by the same transformation. The very important issue, which of the conformal frames, 
Jordan or Einstein, is the physical one, has been the subject of longstanding debate (see [HI W2\ . 
and references therein). In this respect, knowledge of the transitions between these frames is a very 
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useful tool for the construction of new exact solutions. We prove this statement in Sect. [3 

The nonlocal model we will consider is usually studied in the Jordan frame, but recently its 
behavior in the Einstein frame, for the model without matter [311 ES], has been explored too. In 
this paper, we will focus our effort on the study of cosmological solutions of this model, both in the 
Jordan and in the Einstein frames, including the case with matter for the last one. We will consider 
gravity models with a cosmological constant A and including a perfect fluid, and study in detail their 
cosmological solutions with a power-law cosmic scalar factor: a oc t". The solutions thus obtained 
will be proven to generalize some solutions found in [271 (Ml |35]. In the Jordan frame we do not only 
find power- law solutions, but we also prove that other power-law solutions are not exact. We will 
analyze with care the correspondence existing between the solutions obtained in the different frames 
and will demonstrate explicitly how the knowledge of power-law solutions in the Jordan frame can 
be used in order to get power-law solutions in the Einstein frame. 

The paper is organized as follows. In Sect. 121 we start from the action for a general class of nonlocal 
gravity, without specifying the form of the function f{ip), and derive the equations of motion for 
a spatially fiat cosmology in the Jordan frame. The theory is specialized in Sect. El to a model 
characterized by the function f{ip) = /oc"^, and the corresponding power-law solutions are thereby 
obtained. In Sect. Hlwe consider a conformal transformation from the original (Jordan) to the Einstein 
frame, and derive the corresponding equations of motion. Vacuum power-law solutions in the Einstein 
frame are derived in Sect. [51 In Sect. [6l we include the matter sector and obtain the corresponding 
power-law solutions by directly solving the EOM. In Sect. [7] we use the correspondence between 
power-law solutions in the Jordan and in the Einstein frames in order to get brand new solutions in 
the Einstein frame, of a class that cannot be found through the EOM themselves. Finally, Sect. |8]is 
devoted to conclusions. 



2 Nonlocal gravitational action and the equations of motion 

We start by considering a class of nonlocal gravities, whose action is given by 

S = I d^xv^ I + 2A] + , (1) 

where = SirG = Svr/Mpi^, the Planck mass being Mpi = G^^/^ ^ ^ 2 x lO^*^ GeV, while / is 
a differentiable function, which characterizes the nature of nonlocality, with D^^ being the inverse 
of the d'Alembertian operator, A the cosmological constant, and the matter Lagrangian. For 
definiteness, we assume that matter is a perfect fluid. We use the signature (— , +, +, +), g being the 
determinant of the metric tensor, g^^. Recall the covariant d'Alembertian for a scalar field, which 
reads 

□ ^ V^V^ = -La^ iV^g'^'d,) , 

where is the covariant derivative. 

For practical uses, by introducing two scalar fields ip = 0~^R and a Lagrangian multiplier ^ 
(see [2Z]), action (JT]) can be recast as a local action, namely 

Si = y"ci^a;v^|^[i?(l + /(^)) + e(i?-nV^)-2A]+An} . (2) 

Therefore, the original action can actually be regarded as a local one ([2]) in the Jordan frame. By 
varying this action with respect to C, and tp, one respectively obtains the field equations 

= R, (3) 

□e = f,4i^)R, (4) 
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where f,^{ip) = df/dip. The corresponding Einstein equations are obtained by variation of the 
action ([2]) with respect to the metric tensor (yf^,^, as follows 

^g^, [m + d.^d'^ij - 2(A + m)] - R^,m - i {d^id.ip + d^i^d^o + v^9,^ = -k%^)^, , (5) 

where = 1 + f{Tp) and T^^^^y is the energy-momentum tensor of the matter sector, defined as 

_ 2 5{^C,^) 

We note that the system of equations here considered does not include the function ip itself, but 
instead fiip) and f^^piip), together with time derivatives of ip. Also, fiip) can only be determined up 
to a constant, since one may indeed add a constant to fiip) and subtract the same constant from ^ 
without changing the original equations at all. 

In this paper, we assume a spatially fiat Friedmann-Lemaitre-Robertson-Walker (FLRW) uni- 
verse, with the space-time interval 

ds"^ = -dt^ + a^{t)5ijdx'dx^ , (7) 

and consider the case where the scalar fields il){t) and S,{t) are only functions of the cosmological 
time. Thus, the system of Eqs. ([3])-([5]) reduces to 

?>H^^ = - i^^-3/f^ + A + fi:Vm, (8) 

(2H + 3H^^ ^ = ^iij-i/ - 2Hif + A - k^P^ , (9) 

^ = - 3HiP - Q (^H + 2H^^ , (10) 

e = -3Hi-6(^H + 2H'^l4i;), (11) 

where a dot means differentiation with respect to time, t, in the Jordan frame: A{t) = dA(t)/dt, 
and H = d/a is the Hubble parameter. For a perfect matter fluid, we have T(in)oo = Pm and 
T{m)ij = Prndij- The continuity equation is 

p^= -3H{P^ + p^). (12) 
It is useful to add up (IHl) and ([9]), and get 

^ + 5Hij + (2H + 6H^^ ^ - 2A + k^{P^ - pj = 0. (13) 
Note that Eq. (fT3!) is a second-order linear differential equation for 

3 Power-law solutions of the model with /(^) an exponential 
function 

3.1 The model with /(V^) being an exponential function 

Following [3ll |35l [36] , we consider the case where f{i/j) is an exponential function, namely 

= /oe"^ , (14) 
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with /o and a nonzero real parameters. The motivation for considering an exponential function /(V') 
is not only because it is the simplest model with power-law and de Sitter solution^, but also, because 
it is the better studied case among all possible functions /(V') [211 |28l [291 ESI EH ESI ES] (de Sitter 
solutions for this model were discussed in [271 EH ES], and expanding universe solutions with the 
Hubble parameter H = n/t, where n is a nonzero constant, were obtained in [23 ES])- In the present 
paper we will investigate this last type of solutions in detail. 

As advanced in the introduction, we consider matter with the EoS parameter = Pm/Pm being 
a constant but not equal to —1. For power-law solutions H = n/t, Eq. f lT2|) has the following general 
solution: 

Pm(t) = P0t-'"^"-"+'\ (15) 

where po is an arbitrary constant. 



3.2 Solutions with H = n/t 

The goal of this Section is to find the whole set of power-law solutions for the model, characterized 
by the function / given in f|T^ . In this subsection, we present some power-law solutions and, in the 
next one, we will show that no other power-law solutions exist. 

Inserting H = n/t into Eq. flTU]) . the following solution ipit) is obtained, 

a-3n 6n(2n - I) f t 



m = ^it^'^" - 3^_^ ' In J , (16) 

where ipi and to are integration constants. Note that this solution is valid only for n 7^ 1/3 and 
n 7^ 1/2. Consequently, in this subsection the cases n = 1/2 and = 1/3 will be excluded from our 
analysis. We also, specify ipi = 0, so that the function /{ip) takes the following form 

/«'W)^/o(i)"'. ™.-6.!|;^. (17) 

We will show in the next subsection that there is no solution for ipi ^ 0. The cases n = 1/2 and 
n = 1/3 will be considered in the last subsection of the present Section. 

Inserting formulae flTB]) and flTTl) into f ill I) one obtains the following expression for ^ (t) 



, f^V''^ (3^-l)/0 r . 

6 - In , for m = 1 - 3n, 



where ^q, ^1, ^2, and ti are integration constants. 

Furthermore, substituting the solutions described by formulae ( fT5l) . (fT6|) . and (fTSi) . into Eqs. (IE]) 
and (ini), constraints on these integration constants can be obtained 

• For m 7^ 1 — 3n, which equivalently implies the following constraint on the power index n 



3(.-l)±v/3a(3a-2)^ (19) 

3(4a; — 3) 

in this case, we get constraints on the integration constants for A = and A 7^ separately: 



^In models with such solutions, the function f{ip) is either an exponential function or a sum of exponential func- 
tions 1551. 
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— For A = 0, by inserting the solutions (ITSl) . (fT6l) . and (fTS!) into the system (!8|)- (fTT]) . the 
corresponding integration constants are determined to be 



eo = -l 

_ 6{3n — 1 + 3a — 6nQ;)/on^ 6n(2n-l)a/{3n-l) 

(20) 



with to and to be determined by initial conditions, while the power index n is constrained 
by 

2 2a(2n-l) , , 

from which n is expressed in terms of the EoS Wm and parameter a: 



3iym — 6a + 9 ± y (3wm — 6a + 1)^ + 8 (1 — 3«;m) 
6(3wm - 4a + 3) ■ 

This furthermore yields corresponding constraints on the parameters a and w^a for a real 
number n±. 

(3wm - 6a + 1) V 8 (1 - 3wm) ^ . (23) 

Interestingly enough, from Eq. f l22p one finds that, for a ^ 1, one of the power indices 
behaves as n_ — > 1/2, which implies that the Universe asymptotically behaves as a 
radiation-dominated universe, and this regardless of the details of the EoS for the matter 
sector. Thus, we manage to obtain an expanding universe without introducing a cosmo- 
logical constant A. This is because the non-local term D"^/? plays partially the role of a 
dark energy, though it is a decelerating expansion. 

Furthermore, we pay attention to two special values of the EoS parameter: 

(i) When Wm = —1 matter is nothing but just an effective cosmological constant. There- 
fore, this case corresponds to A 7^ 0. 

(ii) If the only matter is radiation, namely, Wm = 1/3, Eq. (!22|) leads to n_ = 1/2, which 
should be excluded. Therefore, the power index is n = n+ for the radiation component. 

For A 7^ 0, the corresponding integration constants and parameters are constrained as 



m = 2 



6n{n + l)/o 



A ' (24) 



Po 
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where the integration constants ^0 and ^1 are to be fixed by the initial conditions, while 
the power-index n is determined by w^'- 

n = — ^ . (25) 

Here we note that, by recalling the definition of m in Eq. (fT7|) and using Eq. (l25l) . one 
finds that the parameter a is constrained by m = 2: 

a=f- <\ . (26) 

2(3wm - 1) 
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Thus, unhke the situation for vacuum solutions with a nonzero cosmological constant 
in [35], here we find that with a matter sector the parameter of the model a is fixed by the 
EoS of matter in Eq. (!26|) . In this sense, the model is spoiled since, in this case, it cannot 
yield a smooth evolution of the Universe for different stages with a given parameter a. 

From fl26|) it follows that there is no power-law solution for non-vanishing cosmological 
constant with matter whose EoS is Wm = 1/3 or Wm = 1, if m 7^ 1 — 3n. The case 
Wm = — 1 corresponds to the cosmological constant as the matter part. Hence, to obtain 
the solutions in this case we should put p(t) = 0, which corresponds to ^0 = —1 from 
Eq. (El]). In this case, from the constraint m = 2, we also obtain two branches for the 
power-index: 

_ 3(a - 1) ± V9a2 + Qa + 9 

~ 12^ ' ^ " 

which coincide with Eq. (25) of Ref. [35], as expected. So the solutions obtained here con- 
tain the vacuum case where the universe asymptotically behaves like a radiation-dominated 
one for large a. 



For m = 1 — 3n, the power-index n is determinecil by the parameter 



a: 



3(a - 1) ± J3a(3a - 2) , , 

For A = 0, one finds power-law solutions with the following constraints on the integration 
constants 

^2 = — 1 , ^0; ^1 are arbitrary parameters, 

_ 3n{n-l)fo (29) 

Po h 



K 

while, again, n is determined by Wm- 



n = 7— ; (30) 

3Wm 



thus, similarly as in the former case for m ^ l — 3n with nonzero cosmological constant A, here 
the parameter a is also fixed by w^: 

^^3(1-^ 
2(2 - 3wm) ' 

which again implies that in this case the model cannot yield the different stages of the Universe's 
evolution with a fixed a. 

For some special values of the parameter a there exist additional solutions: 

— For a = 2/3 and A = 0, there exists a solution with integral constants ^2, ^05 ^.nd ti to be 
specified by the initial conditions, while n, Wm and po are fixed as follows: 



3(6 + 1) 

3 ' 



1 ^ "^^^2 -r -l; 

n = l, w,n = - - , Po = Z2 ■ (22) 



^The equation is the following: 



(3n-l)^ 
Q: — = (J. 

6n(2n — 1) 
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— When a = 6/5 and A 7^ we obtain a solution with ti undetermined while 



^=-^, Po = ^^ 6=-l- (33) 

Note that for a = 6/5 we also have a solution (l29ll with A = and n = 5/9. 

Thus, we have obtained corresponding solutions for both nonzero and zero values of A. They 
generalize the ones previously found in the absence of matter [35]. All solutions for the case m = l—3n 
are new. 



3.3 Proof of the absence of power-law solutions in the case ijji 

Since the goal of this paper is to find all solutions that correspond to H = n/t, in this subsection, we 
consider the case ^/^i 7^ to either obtain new solutions or else rigorously prove that such solutions 
do not exist. 

In the case ipi 7^ 0, Eq. (ITTl) has no solution in terms of elementary functions, thus it is more 
convenient to solve Eq. get ^'(t), and substitute ^{t) = -^{t) - f{i^{t)) into Eq. ([11]), to then 
check if there exist values of the parameters for which the obtained ^(t) satisfies Eq. fill I) , or not. 
The type of solutions of Eq. f[T3]) depends on the value of n, thus we consider some different cases. 

For n ^ —1 and n 7^ —1/3, Eq. f[T3]) has the following general solution 

V|/(t) = + C2t^-3" + 7 ^ rt' - /° ^^"^ (34) 

[n + l){3n + 1) {3nWra — l){n + SnWm — 2) 

where Ci and C2 are integral constants and w^a is chosen so that {3nWui — l)(n + SnWm — 2) be not 
equal to zero. 

In the case n ^ —1 and n 7^ —1/3, the solutions of Eqs. f[TO]) and f[T3l) are given by formulae flMl) 
and ( fT6l) . respectively. Substituting the ^(t) = "^(t) — /{ipit)) obtained into ( fTTl) . we realize that this 
equation is not satisfied. In particular, the non-matching expression is proportional to 



exp 



■/3 



where /3 and 7 are combinations of constants. The exponential term disappears only for ipi = 0. 
Therefore, there is no solution for ipi ^ and n 7^ —1, n 7^ —1/3. Similar calculations show that, in 
the cases n = —1 and n = —1/3, solutions are absent, as well. 

We thus have found all solutions which correspond to H = n/t. Note that, in contradistinction 
to the papers [271 ES], we here include an ideal perfect fluid in the action ([2]) and do not impose any 
restrictions on the parameters and integration constants. 

3.4 Special values of the power index n 

Let us consider the case n = 1/2, which corresponds to i? = 0. Solving Eqs. fITO]) and f lTT]) . we get 

m = i^^t-'^' + ^4, m = ^t^'/' + ^4, (35) 

where ips, ip^, ^3, and ,^4 are integral constants. Straightforward substitution of these functions and 
H = 1/ (2t) into Eqs. ([8]) and ([9]) gives solution for A = with the following conditions on the 
constants: 

^3 = 0, ^4 = -1 - /oe"^* + ^«:Vo, ^- = ^' (36) 
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while po; and ^3 are to be determined by the initial conditions. 
When n = 1/3, Eq. (1T0|) has the following solution: 



(37) 



where tp^ and ^2 are integration constants. The function ^ (t) as a solution of Eq. flTT]) can be given 
in terms of quadratures only. At the same time, solving Eq. ( I13p . we get in terms of elementary 
functions. Thus, if a solution exist, then the corresponding ^{t) should be an elementary function 
as well. We therefore arrive to a contradiction, what proves the absence of power-law solutions with 
n = 1/3. 

3.5 Brief summary of the solutions in the Jordan frame 

To make it easier for readers to look at the whole set of solutions, we list all the solutions we have 
found in this section in Tables [1] and [21 We note that in the case n ^ 1/2, 1/3, the solution for ilj{t) 
is uniquely given by Eq. flTB]) with ipi = 0, i.e. 



6n(2n 



3n 



1 V^o 



(38) 



and thus, this expression is not repeated in the table. It should be noted that for n = 1/2 the 
solution is given by Eqs. fl35|) and fl36|) . while no power-law solution exists for n = 1/3. 



Table 1: Solutions in the Jordan frame for m 7^ 1 — 3n 



m 7^ 1 — 3n 



solutions 



constraints 



A = 



l-3n 



Pm{t) 



1 + 



(3n-l)/o /t 
3?2 + m — 1 yto 

3a(l-2n " 



3n- 1 



6a.n(2n — l) 
^ 371-1 ^ — 3n(li!m + l) 



Eq.(l22D 
Eq.dlHD 



A ^ 



to J 3n + l \to^ 



Pm{t) 



3n\l + Co 



0j^-3n(i«ni+l) 
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Table 2: Solutions in the Jordan frame for m = 1 — 3n 



m = 1 — 3n 



solutions 



constraints 



A = 



at) = -1 + /o(3n - 1) ( ^ ) In ( ^ 



Pra{t) 



t \ / t 



3/o»-(l - ^) ^3n-l_^-3nK„+l) 



Eq. ([28]) 
Eq. ([30D 
Eq. dSI]) 



A ^ 



Pm(^) = 0, 



In 



a = 6/5 
Eq.dSSD 



4 Action and equation of motion in the Einstein frame 

4.1 The Jordan and Einstein frames 

Once a modified gravity theory is recast into its scalar-tensor presentation, it immediately follows 
that both the Jordan frame (where the matter sector minimally couples to gravity) and the Einstein 
one (where the Ricci is linear but matter couples to gravity non-minimally) are available (for a 
description see |1Q]). These two frames are related by conformal transformation 

9,. = ^'gif , (39) 

where we denote the metric in the Jordan frame by g^j^^, while the one in the Einstein frame is labeled 

One soon realizes that the conformal transformation connecting both frames cannot be simply 
interpreted as a coordinate transformation of the theory, and this is the reason why there has been 
a long debate on which of these two frames is 'the physical one', regardless of the fact that the 
mathematical equivalence of the two frames is quite clear [HI |42l |43] . Recent researches have further 
clarified that, at least at the classical level, the two frames are physically equivalent, what means 
that all observational quantities (e.g., the redshift z) should yield the same value in both cases [45j. 
Based on the corresponding solutions in both frames, one can furthermore probe the equivalence of 
frames in the framework of nonlocal gravity inspired models. It is of interest to check the precise 
behaviors of the corresponding solutions in the Einstein frame to see if and, being the case, how much 
they differ from those obtained in the Jordan frame. Moreover, we show that the transitions between 
these frames is a useful tool for the construction of power-law solutions in the Einstein frame. 

4.2 Conformal transformation 

In this Section, we study the nonlocal model ([2]) in the Einstein frame with a perfect fluid. It has 
been shown that those fluids have the same EoS in both frames [20], which implies that the matter 
sector remains unaltered in both. Hence, it will be very interesting to trace the behavior of the 
cosmological solutions in the Einstein frame, with the same matter fields. 

On the other hand, it is known that a theory with higher derivatives in the action will often 
suffer from the ghost problem, i.e. a wrong sign in front of the kinetic terms, which will cause 
an instability problem thus making the theory physically irrelevant. Therefore, it is important to 
examine if the theory contains a ghost or not. To see this non-perturbatively, one needs first make a 
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conformal transformation of the metric to bring the action into the form of the one in the Einstein 
frame [331 EH |35], namely the conformal frame in which the gravitational part of the action ([2]) 
becomes purely Einsteinian. Note that the matter field is assumed to be minimally coupled to 
gravity in the Jordan frame, as given in action ([2]). Power-law solutions in the Einstein frame for 
the model given by the action without matter were considered in [511 ES] • In the following, we 
denote all quantities in the Einstein frame by adding a tag to the corresponding ones in the other 
frame, in order to avoid confusion. 

Let us consider the conformal transformation fl39|) . Using g'^'^ = Q' 

-2gME)^ one obtains the 

relationship between the Ricci scalars in the two frames 

R = - 6 (n^^) Infi + (^^-{^)vf ) Inf^Vf ) Inil)] . (40) 

Inserting this into action (j2]), one immediately identifies the conformal factor as [33 | l34 l [35] 

f^~' = ^ = i + /w + e, (41) 

and then, by introducing a new field 0, 

0^1nl]= -iln(l + /(V')+0= -^InW, (42) 

to remove the Lagrangian multiplier ^, this finally results in the following action in the Einstein 
frame: 

S = l rf^xv^ - 6V-(^)0Vf - 2V-(^)0Vf ^ ^^^^ 

- e2^/,^(^)V^(^VvfV-2e^^A] +e"<^/:^(Q; e^'^g^''^)} . 

In what follows we will derive the corresponding equations of motion in the Einstein frame by varying 
this action. After solving them we will discuss their specific cosmological behaviors. 

4.3 Equations of motion 

By varying action fll02p with respect to the metric obtains the corresponding Einstein 

equations, as follows 

Rlf - IglFJ [R^""^ - 9^^^''^ {Gd.<pdp<P + 2d.<pdp^ + e'^Ud.^dpi^) - 2e^<^A] - ^^^^ 
where we recall that is defined by ( H2l) , and the energy-momentum tensor in the Einstein frame as 



2 5(fi^V^=^£n., 

The field equations read 

□ (^)(60 + V;)-e2%(7'^'^(^)9^^a,z^-4e^^A + «:V^(^4£^ + ^^ = 0, (46) 
2n(^V + 2e2^/,^n(^V + ^7'^''^''^e2<^(4/,v,9^V"9.0 + /,^^a^V"9.V') = 0. (47) 
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At first sight, one may guess that the last term of Eq. P6l) could be troublesome. However, this 
term can be substituted by a combination of the conformal factor and the trace part of the energy- 
momentum tensor in the Einstein frame. To achieve this, from Eqs. ( !39|) and (H5|) . we get 



dg'^'^ d(j) dg 



-2fi-2^'^'^(^)^ = -2^^^^^, (49) 



while from the definition of the energy-momentum tensor in the Jordan frame, we obtain 

r; = - = 4£„ - . (50) 

Inserting Eqs. ( HHj) and (H9|) into (150|) . one recovers the last term of Eq. (H6l) under the form 

AC^ + ^ = . (51) 

4.4 The FLRW metric 



As is well known [20], conformally fiat metrics are mapped into each other. The FLRW metric is 
conformally fiat, so starting from the FLRW metric in the Jordan frame we obtain the corresponding 
FLRW metric in the Einstein one. This leads us to directly start from an FLRW metric with cosmic 
time in the Einstein frame 

ds'^ = - dt% + a%{tE)Sijdx'dx^ , (52) 

where in dtE and aE{tE), the index E denotes the corresponding quantities in the Einstein frame. 
We get 

dtE = n-'dt, aE = n~'a, HE = ^^^^ = n(H-^] (53) 

dtE \ i ' / 

Under the conformal transformation (!39l) , the energy-momentum tensor of a perfect fluid trans- 
forms as [m ED] 

T,^(^) = diag(-p£;, Pe, Pe, Pe) = diag(-p, P, P, P) . (54) 



Using this equation together with Eq. (|53l) . we obtain that the continuity (conservation) equation ( IT2l) 
transforms into the following one in the Einstein frame 

p'j, + ?.He{pe + Pe) = ^(Pe- ^Pe) , (55) 

where the prime denotes derivative with respect to the cosmological time in the Einstein frame, i.e. 
' = d/dtE- 

From ( 15^ . we obtain the EoS parameter Wm = Pm/Pm = Pe/pe, therefore, the conservation law 
(155|) can be rewritten as 

p'e + Pe [SHEil + w^) + 0'(3w;^ - 1)] = . (56) 

It should be noted that if matter just reduces to radiation {w^ = 1/3), then the conservation laws 
in the Jordan and Einstein frames coincide. 
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For the model with a perfect fluid, Eqs. (j4 



k'2 



B7|) acquire the following form, in the FLRW metric, 

2<j> 



L'2 



- -yM" - e''^ 



K^PEi'iWm - 1) 
KV£;(3Wm - 1) 



(57) 

(58) 
(59) 
(60) 



where the first two are independent Einsteinian equations, while the other two are scalar field equa- 
tions. Thus, the complete set of equations is given by ( I57|) -( l60i) . Eq. ( l56i) follows from this system. 
Combining Eqs. ( 157|) and (l58l) . one obtains the remarkably useful result that the equation for a scalar 
can be reduced to an algebraic equation, as follows. 



Ae4<^ = H'+ ml + 



k'^Pe 



Wn, - 1' 



For A 7^ 0, one can formally obtain from fl6Tl) the expression for the conformal factor ^(t^;) 



0(t, 



In 



(61) 



(62) 



It should be also noted that, by combining Eqs. (158!) and (!59|) . one can eliminate the fiip) and matter 
terms to obtain a second-order differential equation in terms of ip{tE) and (f){tE), namely 

60" + tp" + 3He{64>' + tp') + 6{H'e + 2Hl) + 60'' + 20V = . (63) 

It is most convenient to derive the expression for ipltE) from this equation once He and 0(t£;) are 
known. 



5 Power-law solutions in the Einstein frame for the model 
without matter 

We here investigate power-law solutions (with He = UE/tE) in the Einstein frame. We consider the 
case of an exponential /(V'), given by ( IT^ . For the model with pe = 0, the cases with and without 
matter will turn out to be essentially different. 

If A = 0, from fl6T|) we obtain the following power-law solution 

H'e = -3HI He= ^^^^\^^y (64) 

meaning that all solutions correspond to the power-law behavior of the Hubble parameter with 
He = 1/3. Thus, the general solution of the system fl57l)- fl60l) . which corresponds to arbitrary initial 
conditions for the scalar fields, yields just one cosmological evolution, specified by flMl) . 
In the case where A 7^ and matter is absent, from flB21) we obtain 

(65) 

Note that there is no solution with n£; = l/3ifA7^0. Eq. ( l63l) is a linear differential equation for 
ipitE)- Inserting the function 4>(tE) we just obtained into this equation, we get the following solution 

i^itE) = V^o4"'"" + rriE In (j-^ , (66) 



0(t, 



I- 



ue^Sue - 1) 
Atl 
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where niE = 3(2ri£; — l)(4n£; — 3)/[2(2 — Sn^;)], while ipo and to are two integration constants. Inserting 
fl65|) and (166|1 into the system (1571) - (1601) . one gets the following constraints on the integration constants 



aniE 



0, 



1 



to = 3/o(2r2s - 1) 



He 



K'iriE - 1) 



a + 



2{'inE - 2) 







(67) 



(6J 



3(2ns-l)(4ns-3) 

from which the two branches corresponding to the index ue are expressed in terms of the parameter 
a, as 

_ 15a - 3 + v/3(3a2 + 2a + 3) 

' (69) 

15a - 3 - ^/3(3a2 + 2a + 3) 



n 



E{1) 



24a 

We note that this solution is the same as the one found in [35] , where the vacuum solution with non- 
vanishing cosmological constant A is constructed by conformal transformation from the corresponding 
one in the Jordan frame. Also we note that for any range for the parameter a, the range for nE{i) is 
nE(i) e (1/2, 3/4), while for nE{2) we have nE{2) G (— cxd, 1/2) U (3/4, + oo). 

There are a few special cases of the parameter ue for which the above mentioned solution does 
not exist. In these cases the parameter m^; is either equal to zero or it does not exist. Let us consider 
all of them in detail. 

• If n^; = 1/2, then 

^(ti^) = Ci + C2V^. (70) 

Substituting this solution into Eq. (!57|) . we conclude that this equation is not satisfied whatever 
be the constants ci and C2. 



If He = 2/3, then 



12 Vt 



+ Ci In 



tE 

f 



(71) 



Substituting this solution into Eq. we conclude that this equation is not satisfied for any 
value of the constants Ci and T. 



• If n^; = 3/4, then 

^{tE) = Ci + C2t-E^\ (72) 

Substituting this solution into Eq. ( 157|) . with the condition C2 = 0, from Eq. ( I60l) we conclude 
that this equation is not satisfied for any value of ci. 

In summary, we have obtained power-law solutions in the Einstein frame for the model without 
matter, as follows, 

Ue 



H{tE) 
<P{tE) 



4 



a 



UEjSnE - 1) 
Atl 

tE 

3fo{2nE - 



(73) 
(74) 



lAjSuE - 1) 
Ue 



(75) 



where the parameter ue is connected with a by fl69|) . Note that ue ^ 1/2, ue 7^ 2/3, ue 7^ 3/4, 
and Ue 7^ 1/3. This is valid for A 7^ 0. If A = 0, then all solutions of this system have power-law 
behavior with n^; = 1/3 and note that the solution obtained has no arbitrary integration constant. 
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6 Power-law solutions in the Einstein frame for the model 
with matter 



6.1 The case A = 

For H = UE/tE, from Eq. fl^ one gets that 

PE{tE) 



2nE{3nE - 1) 
- w^)t% 



Inserting this into Eq. ( l56l) . the expression for (f){tE) follows 



(76) 



(77) 



where ti is an integration constant. Formula f l77|) is vahd for ue 7^ 1/3, 7^ 1/3, and 7^ 1. 
Furthermore, using Eq. (!63|) we obtain the solution for ip{tE)'- 



^1 



tE\ ^ ^ l2{2nE - 1) [3(1 - w^) + UEi'iw^ - 5)] /^t^ 



U 



[5 - 3Wm + ME^Wra - 3)] (3Wm - 1) 



(7J 



where ^2 and ipi are two integration constants. Inserting Eqs. fl76|) . fl77|) . and flTHl) into the sys- 
tem (I57|) - fl60|l . one obtains the following constraints on the integration constants and parameters 



nE± 

h 



0, 



3 [(3w^ - 11) + 2a(ll - 9w^)] ± |1 - 3wJ^/3 [3{2a -w^ + 1)2 - 16a] 



6 [3{w^ - 1)2 + 4a(5 - 3^^) - 12] 
{3nE - l)(3wm - 1) 



6(2ni5-l)(w7„,-l)/o 



Sium — 1 



(79) 
(80) 



(81) 



It follows that Eq. (IHOj) gives rise to constraints on the parameters a and for nE± to be a real 
number: 



3(2a - + 1)^ - 16a ^ 



(82) 



Interestingly enough, in the limit a 1, using Eq. f lHOj) . we obtain the following asymptotic 
behavior for the index nE± 



nE±{a > 1) 



11 — Qwyn ± 11 — 3Wt, 



3) 



4(5 - 3w 

One immediately realizes that, for = —1, nE+ — > 0.75 asymptotically, what corresponds to an 
upper bound on the value for the index in the vacuum case, in the Einstein frame (see Eq. (38) of 
Ref. p5]). In fact, when w^a = —1, the solutions (ff6l)-(j8T]) reduce to the following form: 

nE{3nE - 1) 



pE{ti 



0(ts) = -^ln(^) , V^(t^) = -ln 



3nE — 1 
3foti{2nE-l) 



4) 



while the corresponding energy density in the Jordan frame is just the cosmological constant A. By 
using conformal transformation Eq. fl5^ and comparing it with the solution for pE{tE) in fl8^ . we 
can express the integration constant ti in terms of A, as 



nE{3nE - 1) 
fi:2A 
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^5) 



It is straightforward to see that Eq. (!84|) recovers the vacuum solutions, Eqs. ( I7i|) - (!75|) . 

Another interesting asymptotic behavior is the one for Wm > 1/3, nE+ — > 1/2, while for < 
1/3 we get He- — > 1/2, regardless of the value of the EoS parameter w^a- 

We thus have derived solutions for all nonzero values of n^;, but for he = 1/3. Also, we assume 
that 7^ 1/3 and ^ 1. Solutions in these special cases will be considered in the next section, 
by conformally transforming the power-law solutions obtained in the Jordan frame into the Einstein 
one. 



6.2 The case A 7^ 



6.2.1 



In this case, we cannot obtain a general solution for the system flFr|) - flBn]) . Nevertheless, from Eq. fl521) 
one finds that = 1 is a special case where the conformal factor ^(t^;) has already been solved. 
Thus, in the following we shall consider the case = 1. 

From Eq. (162|) one obtains the expression of the conformal factor: 



<t'{ti 



4 



{Sue - l)nE 



Atl 



and the condition ue 7^ 1/3. Substituting this function into Eq. ( 156|) . we get 



PE{tE) oc if: 



l—GriE 



On the other hand, inserting the expression fISB]) into Eq. 

3(2ns-l)(4nB-3) 



2(3ns - 2) 



In 



tE 



ip(tE) can be solved as 



^ ,2-3nE 



^6) 



^7) 



(88) 



where Ci and Ti are arbitrary constants and nE 7^ 2/3, ue 7^ 1/2, ue 7^ 3/4. 

However, substituting this function into fl57j) . we discover that this equation cannot hold for any 
value of Ci and Ti. Note that for ue = 1/2, ue = 2/3, and ue = 3/4, Eq. fl63l) has different 
solution^, but the system of equations fl57|) - fl60|) is not satisfied for these values of rtg, either. Thus, 
we conclude that there is no consistent solution of Wm = 1. 



6.2.2 Wrny^l 

If 7^ 1, then we can express via pE{tE) using Eq. fl62l) and, furthermore, obtain the following 
master equation for pEitE) from Eq. fl56l) : 

Mpe) + A2{pe) ^ Q .gg. 

AtE[K.^PEtl{w^-l) + 2nE{?>nE-l)] ' ^ ' 

where, to simplify the expression, we have defined Ai{pe) and A2{pe) as the following two functions 

of PE{tE) 

Ai{pe) = AhePe [3k^(«;^ - l)pst| + {Me - 1) (6n£;(l + w^n) + 1 - 3^^)] , (90) 
A^ipE) = p'EtE[M\w1-l)pEtl + ME{'inE-l)\ . (91) 

For generic values of the parameters and he the general solution cannot be found in terms of 
elementary functions and it is only possible to cast this equation in transcendental form. At the same 
time, it is actually easy to solve it for some particular values of the constants, as the ones which 
follow. 

^These solutions coincide with d7ni)" ((7^ . 
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For ue = 1/3, we obtain from Eq. fl62l) : 

while the master equation Eq. (189|) has the general solution 

= Pat'E^^- (93) 

A straightforward calculation shows that the system (l571) - (!60ll has no solution for ue = 1/3 
if Wm 7^ ±1. (Note that we always assume that ^ —1, because matter with = — 1 
coincides with the cosmological constant.) 

riE = 1/2 



In this case, solving the system ( 157|) -( 160|1 . it is found that there only exists the trivial so- 
lution ipEitE) = constant. Recalling the action ([1]) we are considering, such a constant solution 
just corresponds to a rescaling of the Ricci scalar R, or equivalently, to a rescaling of both the 
Newtonian constant G = Svt/k^ and the cosmological constant A. 

The correspondence of the special solutions considered above is that in Einstein frame, the case 
of a nonvanishing cosmological constant dose not admit a power-law style matter sector other than 
Pe oc t^^ for the scaling solution He oc nE/tE- To see this clearly, we consider the ansatz pEitE) = 
Pit^E- Then the master equation for the matter sector, Eq. fl89|) . reduces to the following form 

AnE{3nE - 1) [1 + 2p - 3w^ + 6nj;(l + w^^)] + SpmH^j^^jAnE + p){wi - 1) ^ ^ 
4tE[K^Pit''E^\w^-l) + 2nE{3nE-l)] 

Let us first consider the case p ^ —2. Eq. (1941) is then satisfied if and only if each of the two 
terms in the numerator are equal to zero. For the first term to vanish it is required that0 He = 1/3 
or 1 + 2p — 3w^ + 6n£'(l + Wm) = 0, while for the second, that p = —^ue (recall that = ±1 is 
excluded here), from where the only allowed possibility here isn£; = l/3, p = —4/3. As discussed 
above, this does not satisfy the EOM flSTjl - fl^ . 

Solutions with p = —2 correspond to the ones already considered in Sect. 17. 4[ 



7 Relationship between power-law solutions in the Jordan 
and in the Einstein frames 

7.1 Conformal transformation between power- law solutions 
7.1.1 General expression for the conformal factor Q 

In Sects. [5] and [6] we have found solutions for matter and massless sectors respectively. However, 
we recall that in Sect. 16. H general solutions were found except for some singular values of n^; and 
Wm, since it was not possible to separate variables in the system (I57|) -fl60l) to explicitly solve for the 
functions ip and (p in these special cases. 

Meanwhile, in the vacuum case, it has been found that some power-law solutions in the Jordan 
frame correspond to other power-law ones in the Einstein frame, by using conformal transforma- 
tion [35]. Here we will generalize this correspondence for the model with a matter sector, and 

"'We always assume that ue ^ 0. 
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furthermore use the conformal transformation to obtain those special solutions in the Einstein frame 
which are very difficult to obtain by directly solving the system (l57|l - (!60!l . 

First, we formulate the differential equation for the conformal factor under which power-law 
solutions in the Jordan frame correspond to other power-law solutions in the Einstein frame. Using 
(155]) . we have 

he = '^ = - m , (95) 

lE t 



where we recall that ' = d/dt. This immediately gives a relationship between the cosmological times 
in the Jordan frame, t, and in the Einstein one, Ie, as follows 



tE 



riEt 



nfi(t) - n{t)t 



Taking the derivative of the equation above with respective to t and inserting dtE/dt 
obtain the following differential equation for the conformal factor Q{t): 



Q - 



n^t 

This equation has the following general solution: 



n-Et^ 



0. 



(96) 



we 



(97) 



nit) 



{j^Y Bo (f^)" " + Bi 
(i) [i3,ln(|,) + l] — 



l-n 



n = 1, 
n = 1, 



riE 7^ 1, 

riE 7^ 1, 
Ue = 1, 

He = 1, 



(98) 



Bif^ exp 

(*)"■. 

where Bq, Bi and T are integration constant^. 

In principle, as it was done in Ref. [35], one can now obtain power-law solutions, by using Eq. (pHjl . 
to conformally transform the corresponding solutions from the Jordan frame into the Einstein one. 
However, difficulties arise if one directly wants to obtain a general solution as 



7.1.2 Conformal factor corresponding to power-law solutions in the Jordan frame 



As already stated above, although we have obtained the general solution (jHE]), the complication of 
this expression renders the analysis of the solutions rather difficult. Another approach to construct 
the conformal factor Q is to insert solutions found in the Jordan frame into Eq. (1411) to obtain the 
expression directly case by case. Let us consider the case m 7^ 1 — 3n, i.e. use the corresponding 
solutions ffT6D and ffTSD: 



m = ^o + ^i{^ 

^0 



l-3n 



+ 



(3n-l)/o ft 



3n + m — 1 \tQ^ 
Thus, one obtains the expression for fl by employing Eq. fj4T]) 

^]-2(t) = 6-1^ + c2t'-3" + i + eo =^ n{t) 



, / X "^1 ft 

^(t) = _ln 

a \to 



1 



(99) 



(100) 



^Cit"' + C2tl-3" + 1+^0 

where, for simplicity of the expression, we have defined the two parameters Ci and C2 as follows 

C^i - f^^^/oto™ , C,^^,tr-\ (101) 
3n + m — 1 



^Note that we use Bq only to include the case Bq = 0. Nonzero Bq we always can put equal to unit. 
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and again we recall that m is given by f lT7|l . Observing Eq. fllOOp . we realize the constraint on 
the integration constant = since otherwise we can never find a correspondence of power-law 
solutions between the two frames. Recalling now Eqs. ( 120|) and f l2^ . one immediately realizes the 
implication of this constraint: there is no correspondence between power-law solutions in the case of 
non-vanishing cosmological constant A and the matter sectors. 

It is easy to see that there are two simple cases, namely Ci = and C2 = 0, when the power- 
law Hubble parameter in the Jordan frame H, after conformal transformation, yields a power-law 
function in the Einstein frame. We consider these cases separately in the following subsections. In 
Subsects. 7.2-7.4 we will consider the case A = 0, whereas the case of nonzero A will be dealt with 
in Subsect. 7.5. 



7.2 The case Ci = 



The case Ci = is special and corresponds to fixing the index n. Indeed, from Ci = 0, Eq. (llOip 
yields 



n = riQ 



-6 + 3a ± v/3a(3a - 4) 



6(2q; - 3) 

Comparing fl22|) with (11021) . one realizes that = 1. Inserting Eq. (llOOp into (1531) . we obtain 



(102) 



3(1 -no) 



3(l-"o) 



-1 



'Cot— 



'Co 



3(1 -no) 
2^/C'o 



1 — JTlQ 

3(l-no) 



(103) 



hence, using the relationship (153|1 . we find ue = 1/3, thus the Hubble parameter in the Einstein 
frame is 

HE{tE) = ^. (104) 

6tE 

Using the definition = lnf2, we obtain the expression 



3^0 - 1 , (tE 
In 



3(1 -no) \t^ 
Inserting (fTM]) and (fTO into Eq. ([M]), V'(^s) is found to be 



3(1 -no) 



(105) 



2+^ 



I ,'tj;V-«("o-i) , 4no(2no-l) , As 
^^'^^-^^^Z) +(no-l)(3.o-l)^Hr 



(106) 



with two integration constants ifj^ and t*. On the other hand, from ( l99l) and f ll03p . we can rewrite 
the solution (fT6l) in terms of the variables in the Einstein framej^ 



^{tE) = -\n\- 
a \ to 



3(1 - no; 



tE 



3(1-710) 



4no(2no - 1) , ( tE 
m 



(no-l)(3no-l) \l 



(107) 



where we define the constant ti as follows: 



2JC'o 



3 (1 - no. 



3(l-no) 
-t 2 

''0 



(108) 



^ Again we recall that the integration constant ■01 — 0. 
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Comparing Eqs. (11061) and (11071) . one immediately finds that ip^ = and t^, = t4. Thus, we have 
obtained a new solution for he = 1/3 and Wm = 1 in the Einstein frame by using conformal 
transformation, which would be difficult to obtain directly from the system (I57|) - (l60|) . 

To summarize, by using conformal transformation, we have found the following, new solution in 
the Einstein frame: 

1 



HE{tE) 



Me' 

SUq — 1 

3(1 -no) 
4no(2no 



In 



(no - l)(3rao - l] 



In 



tE 



(109) 



with EoS = 1 and a constant constrained by fl571) - fl60l) . as follows 



P3 



8^0 /o (1 - 2n 



2(3no-l) 
3(no-l) 




;iio) 



where and ^4 are constants defined in Eqs. (11051) and (11081) . respectively. Note that no is not a 
free parameter, because it is connected with a. Solutions have been found for arbitrary nonzero a, 
except foi0 a = 4/3, which corresponds to no = 1. 



7.3 The case Co = 



In the case C2 = (or, equivalently, ^1 = 0), similarly as in Sect. 17.21 the relationship between t and 
tE can be obtained as 

Using this relation, the conformal factor can be expressed in terms of the variables in the Einstein 
frame, namely 

2m 

/ »r) _|_ 9 \ m+2 

n-^(tE) = C^lZ^tE] . (112) 



Inserting (11121) into Eq. ( l53l) . we obtain the Hubble parameter in the Einstein frame 

HEitE) = + ■ (113) 
(m + 2)tE 

Thus, by setting the integration constant ^1 = 0, we obtain the correspondence between the power- 
law solutions in the Jordan and Einstein frames, and identify the index in the Einstein frame with 
the corresponding one in the Jordan frame, as follows 

m + 2n [3(2q; — l)n + 1 — 3aln 

riE = = — n ^ —■ (114) 

m + 2 6an2 - 3n(l + a) + 1 ^ ^ 



^The values of no equal to uq ~ 1/3, which corresponds to a = 0, and uq = 1/2, which does not correspond to any 
finite value of a, are also excluded (see Sect. 3). 
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<P{t,) ^lnn= 2-3Ml + ^m) 1^ / 1,. ^^^^^ 



We recall that the parameter n is determined by a and in the constraint fl22!) . For practical use, 
we rewrite this constraint as 

_ (3n - 2 + 3wm^)(3?2 - 1) 
" ~ 6n(2n - 1) ' ^ ' 

and insert it into Eq. flll4p to express n in terms of Ue and lym! 

„ = ^P'"^ - 1) . (116) 

3nE{Wra + 1) - 1 - 3Wm 

Inserting (11160 into Eq. ( 11120 . we can also obtain the expression for the scalar field (f) 

3ng(l + Wm) / 1 

where is an integral constant defined by 

h = T' (118) 

(m + 2)Cf 

hence connected with to and /o by Eq. ( llOip . We see that the expression ( I117P coincides with ( 177|) 
and thus, generally speaking, we get only solutions previously obtained in Sect. 6.1. Nevertheless, for 
some values of the parameters we do not get the solutions of Sect. 6.1. Let us check the possibility 
to get these solutions using the conformal transformation. 

Recall that in Sect. 6.1 we did not get solutions for Wm = 1/3, = 1, and ue = 1/3. If = 1, 
then from Eq. (!6T]) it follows that ue = 1/3 and vice versa. Thus, substituting = 1 into (I117p . 
we get 4>(tE) = 0. Straightforward substitution into Eqs. ( 157|) - (!60|) shows that there is no solution 
in this case. 

For Wra = 1/3 the formula (I117p is not acceptable. We will consider this case in the next 
subsection. 



7.4 Radiation case 

Up to now we have not obtained any power-law solution in the Einstein frame for the case when 
Win = 1/3, i.e. with a radiation sector. In this case, using the conservation law (156 p . we get 

p(t^) (X t-^"- , (119) 



while Eq. ( 16 ip with A = yields 
which implies that 



3?2£;(3n_E_- 1) 



PEitE) = " (120) 

K i I 



ns = -, PEitE) = ^^. (121) 

However, since (f>{tE) couples with ip{tE), we cannot solve these equations directly from the EOM 
( l571) - (!60l) in the Einstein frame. Recall now that, having obtained the solutions (JT6l) and ( JT8l) in 
the Jordan frame we can, in principle conformally transform both into their corresponding forms 
in the Einstein frame. To achieve this goal, we first notice that, using Eq. (11140 . we can find the 
power-index in the Jordan frame n from the corresponding one in the Einstein frame ue = 1/2, a^ 

n = . ■ (122) 

3 1 — a 



^Note that at — 1/2 and = 1/3, we can not use (|116l) to define n, because both numerator, and denominator 
in this formula are equal to zero. 
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Inserting this into Eq. flll2p . we then get 4>(tE) as 

Afj. \ 1 - 3q; , f tE 



6a - 4 \t2 



(123) 



where t2 is an integration constant. And inserting this solution into Eq. f l63|l . one obtains the solution 

for ip{tE): 



tsJ a (3a -2) V^s 



(124) 



with two integration constants, ip2 and ^3. Eqs. ( [Fr|) -( pU]l introduce constraints on these integration 
constants, namely 



^2 


= 0, 






' 1 


h 




Jo 



(125) 
(126) 



Thus, Eqs. f lT2T|) . f lT23|) - ffT26|) supplement the solutions in Sect. O 



7.5 The case A 7^ 



In the case A 7^ solutions in the Jordan frame are described by ( !T6|) and (fT8|) . To get the cor- 
responding power-law solutions in the Einstein frame we need select the case m 7^ 1 — 3n and put 
^0 = —1- From we get that the system does not include matter: po = 0. It is easy to see that, 
for any nonzero n, Ci 7^ for m = 2, so we put ^1 = and consider the case C2 = 0. From (llUUp . 
we obtain 

n{t) -- ^ 



Therefore, 



and 



n + 1 
2tE 



He 



n + 1 



In(VC^t) 



-In 

2 



2 V CitE 



-In 

4 



{3nE - l)nE 



(127) 



(128) 



(129) 



where Eq. f llOip is used in the last step. So we reobtain solution ( l65ll . Using ( 1271) . we reobtain the 
condition fl69|) . Thus, in the case of nonzero A we can use the power- law solutions of the Jordan 
frame to get the corresponding solutions in the Einstein frame, but this way is here not more effective 
than a straightforward search for power- law solutions of the system (I57j) - fl60l) . 



7.6 Brief summary of the solutions in the Einstein frame 

To render it easier for readers to examine the full set of solutions, in Tables |3l H] and |5] we list all 
those corresponding to the Einstein frame. It should be noted that all solutions in these three table 
correspond to the case A = 0. As discussed in Sect. and 17.51 for the non-vanishing cosmological 
constant case in the Einstein frame, only vacuum solution can be found, namely those in Eqs. (169 p . 
dZH) and (175]). 
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Table 3: Solutions in Einstein frame for he 7^ 1/3, Wm 7^ 1, 1/3 

solutions constraints 

^«.) = '?fr''jy7""'+"f"°'-f' ln 1 . ^ Eq. ia 

[5 - 3wm + 3n£;(wm " 3)J (3u;m - 1) V ^2 



Eq. dHUD 
Eq.dHU) 
Eq.([82]) 



2nE{3nE - 1 



Table 4: Solutions in the Einstein frame for he = 1/3, Wm = 1 



solutions 



constraints (Eqs. ffT02|) and ffTTOj) ) 



3no - 1 ^^f^ 



3(1 -no) V^s/' 

4no(2n - 1) /t^ 
K-l)(3no-l) ''U4y ' 



3{no-l) 



P3 



no 



^no /o (1 - 2no) / 



2(3no-l) 
3("0-l) 



3^2 (1 — no) 



/2 

^4 



-6 + 3a ± V3a(3a - 4) 
6(2a - 3) 



Table 5: Solutions in the Einstein frame for ue = 1/2, Wm = 1/3 



solutions 



constraints (Eq. f ll26p ) 



^ ' 6a -4 



PE{tE) 



3a - 1 Yni^— 



a (3a — 2) \t 
3 



fo 



3a 
T 



8 Conclusions 

In General Relativity, power-law solutions of the type H = n/t correspond to models with a perfect 
fluid whose EoS parameter = -Pm/Pm is related to the power-index by = — 1 + 2/(3r;,). It 
is interesting to try to find similar power-law solutions in modified gravity theories, in order to 
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check how much they deviate from those for GR. In this paper, we consider power-law solutions in 
a class of nonlocal gravity models stemming from the widely probed and quite fashionable function 
/(□^^i?) = /oe"° ^ and which include a perfect fluid with constant EoS parameter 

By recasting the original nonlocal action ([T]) into its local presentation we have obtained 
power-law solutions for this model, with and without the cosmological constant and both in the 
Jordan and in the Einstein frames. In the Jordan frame we have reached the remarkable conclusion 
that all power-law solutions could be find (see Sect. |3]), what is a most interesting outcome of this 
paper. In the Einstein frame, we have correspondingly obtained the power-law solutions either by 
directly solving the EOM, or by performing a conformal transformation of the solutions obtained 
in the Jordan frame. For this purpose, we have generalized the correspondence between power-law 
solutions in the Jordan and Einstein frames as obtained in [35], in order to duely include the matter 
sector. By using this powerful, non-trivial tool, we have obtained the solutions when = 1/3 and 
Win = 1, in which cases it was very difficult to obtain the corresponding solutions by directly solving 
the system ( 157|) -( 160|) . Hence, we have shown explicitly how the construction of solutions by using 
conformal transformation between the two frames does proceed, thus proving that the method offers 
a valuable alternative in the search for new solutions. 

In [38], it has been shown that not only models with exponential f{ip) can have power-law and de 
Sitter solutions. It would be interesting to consider power-law and de Sitter solutions in the models 
where /(V') consists of a sum of exponentials. Another direct generalization of the present analysis 
is to include several components of perfect fluids with different constant values of 

In [5^ I5B] . deSitter solutions in nonlocal models were found. It is interesting to check the 
possibilities for the Universe evolution as obtained from these models, from an inflationary de Sitter 
stage to the late power-law Universe and, furthermore, to check for deviations from the standard GR 
case, and its distinction from other modified gravity theories. 

As is widely known, theories with higher derivatives will often suffer from a ghost problem, 
namely a wrong sign in the kinetic term, resulting in a dangerous instability problem. A good aspect 
in making use of the conformal transformation technique between the two frames is to obtain the 
corresponding ghost-free conditions [331 Ell ES]. It will be interesting to test the solutions obtained 
in this paper to find the constraints on the parameters, hence to check for the possibility to obtain a 
more realistic model which can be responsible for the current observed acceleration of the expansion 
of the Universe. This work is now in process. Also, the analysis of the stability of the solutions here 
found will be performed [IB] . 
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